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NOMENCLATURE 


F ,F  ,F 
X y z 


dimension  of  plate  parallel  to  the  x direction,  in.  (cm). 

dimension  of  plate  parallel  to  the  y direction,  in.  (cm^. 

2 

elements  of  anisotropic  diffusivity  matrix,  in  /hr 
(mm^/sec). 

diffusivity  par^illel  to  tl^  fibers  in  a unidirectional 
composite,  in^/hr  (mm^/sec). 

diffusivity  of  resin,  in  /hr  (mm  /sec). 

diffusivity  normal  to  the  fibers  in  a unidirectional 
composite,  in^/hr  (mm^/sec) . 

rate  of  transfer  parallel  to  the  x,  y,  and  z axis, 
respectively  of  diffusing  substance  through  a unit  area 
of  cross-section,  % weight  gain/hr /in^  (%  weight  gain/ 
sec/mm  3). 

plate  thickness,  in.  (cm). 

% weight  gain 


% weight  gain  at  equilibrium 

3 3 

% weight  gain/in  (%  weight  gain/mm  ). 

3 3 

% weight  gain/in  (%  weight  gain/mm  ) at  equilibrium. 


temperature,  °F  (°C). 


time,  hr(8ec). 

fiber  volume  fraction,  %. 

angle  of  fiber  axis  relative  to  the  x axis  in  a uni- 
directional composite,  radians. 
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SECTION  I 
INTRODUCTION 

Moisture  absorption  in  the  presence  of  a humid  environment  is  a basic 
physical  property  associated  with  current  epoxy  resins  utilized  in  high 
performance  composite  materials.  Such  absorption  causes  plasticization  of 
the  resin  to  occur  with  concurrent  swelling  and  lowering  of  the  resin's  glass 
transition  temperature  (temperature  at  which  the  polymer  changes  from  a 
glassy  solid  to  a rubbery  solid).  Because  of  the  lowering  of  the  resin  glass 
transition  temperature,  there  is  considerable  interest  in  determining  the 
mechanical  properties  of  composite  laminates  at  various  temperatures  in  the 
presence  of  moisture.  Such  experimental  evaluation  requires  a knowledge  of 
the  moisture  diffusion  process  in  order  to  perform  accelerated  moisture  con- 
ditioning and  to  determine  moisture  content  and  distribution  during  elevated 
temperature  tests. 

The  mathematics  of  diffusion  as  applied  to  laminated  composites  is 
summarized  in  the  present  report . Theoretical  results  are  compared  to 
experimental  data  for  graphite- epoxy  composites,  and  some  anomalies 
associated  with  laminates  are  discussed. 


SECTION  n 

DIFFUSION  IN  A LAMINA 


Consider  an  off-axis  unidirectional  lamina  as  shown  in  Figure  1.  For 
the  x-y  axis  system  parallel  to  the  sides,  the  diffusion  equation  takes  the 
form  (Reference  1) 


m,.+F  +F  +F  =0 
t x,x  y,y  z,z 


where  a comma  denotes  partial  differentiation.  The  rate  of  transfer  of 
diffusing  substance  is  related  to  the  concentration  gradient  through  Fick's  law 
for  an  anisotropic  media,  which  takes  the  form  (Reference  1) 


‘^11  ^12 


*^12  *^22 


Substituting  equation  (2)  into  equation  (1),  assuming  the  diffusivities  are  con- 
stant, yields  the  results 


m . = d, , m + 2d,  - m 

,t  11  ,xx  12  ,xy 


+ d m + d_-  m 
22  ,yy  33  , zz 


For  the  case  of  a unidirectional  material  containing  square  of  hexagonal 
fiber  packing 


*^33  ^ ^T 


The  in- plane  diffusion  coefficients  in  equation  (2)  can  be  expressed  in  terms  of 
the  diffusivities  parallel  and  transverse  to  the  fibers  through  the  transformation 
relationships  (Reference  1) 


2 2 

d^  j = dj^  cos  © + d^  sin  © 


*^12  ■ (d^- dj^)cos  © sin  © 


2 2 

*^22  ” 6 + cos  © 


SECTION  m 

MICROMECHANICS  OF  DIFFUSION 

Approximate  micromechanics  expressions  for  and  d^  have  been 
derived  by  Springer  and  Tsai  (Reference  2).  For  many  advanced  composites 

of  interest,  moisture  absorption  occurs  only  in  the  polymeric  matrix  material 

4 

with  the  result  (Reference  2) 

= (1  - 2VVj/ir)djj  (6) 


Combining  equations  (5)  and  (6)  yields 


SECTION  IV 

DIFFUSION  IN  A LAMINATE 


I 

i 

i 


r 


Since  laminates  are  constructed  of  plies  having  fiber  orientations  in 
various  directions,  the  diffusivity  in  any  one  direction  will  vary  through  the 
thickness.  A solution  for  a laminate  would  require  the  solution  of  equation 
(3)  ply- by- ply.  A simplified  approach  can  be  developed  by  using  "effective” 
constant  diffusiviti  s in  conjunction  with  equation  (3)  which  account  for  differ- 
ences in  diffusivities  for  individual  plies.  This  approach  is  analogous  to  the 
development  of  laminated  plate  theory. 

1.  EFFECTIVE  DIFFUSIVITY  THROUGH-THE-THICKNESS 

Consider  a laminated  slab  subjected  to  one  dimensional  diffusion  through- 
the- thickness,  i.e. 

m = m(z)  (8) 

with  the  surface  conditions 

m(-  h/2)  = m^  , m(h/2)  = m^  (9) 


where  m^  and  m^  are  constants.  For  steady  state,  equation  (1)  reduces  to 


z.  z 


= 0 


(10) 


Integration  of  equation  (10)  yields 


F = k,  = constant 
z 1 


(11) 


Equation  (2)  in  conjunction  with  equations  (8)  and  (11)  yield 

F^  = F^  = 0 (12) 

Separation  of  variables  in  equation  (13),  followed  by  integration  through-the- 
thickness,  and  application  of  the  surface  conditions,  equation  (9),  leads  to  the 
relationship 


5 


I 


1 


4 


I 


kj  D = 


where 


/ 


d33(.) 


Effective  diffusivity,  denoted  by  a bar,  is  defined  by  the  relationship 

= ^33 

Integration  of  equation  (16)  in  conjunction  with  equations  (9)  and  (11)  yields 


= ”^1  - "^2 


Combining  equations  (14)  and  (17)  yields 


33  D 


For  a laminate  of  N layers  equation  (15)  takes  the  form 

N , i 


° = E ^ 

• 1 j * 


i=l  d. 


where  the  superscript  i denotes  properties  of  the  ith  layer, 


2.  INPLANE  EFFECTIVE  DIFFUSIVITIES 

Consider  a laminated  plate  of  thickness  h having  infinite  length  in  the  y 
direction  and  subjected  to  the  surface  conditions 

= m,  + (m_  - m,)^  , 


m(x,  ±h/2) 


(20) 


Due  to  the  variation  of  the  diffusivitiea  with  respect  to  z,  the  problem  is  two 
dimensional  in  nature.  A cursory  examination  of  equations  (1)  and  (2)  reveals, 
however,  that  the  steady  state  solution  which  satisfies  the  boundary  conditions, 
equations  (20)  and  (21),  is  the  surface  function  defined  by  equation  (21).  Sub- 
stituting equation  (21)  into  equation  (2)  yields 


(m^  - nij) 


(m^  - m^) 


F = 0 

z 


(22) 

(23) 

(24) 


Integrating  equations  (22)  and  (23)  with  respect  to  z and  dividing  by  the  thick- 
ness, h,  yields 

h/2 


u 

-h/2 

if 

-h/2 


_ (m  - m ) 

F d = d,,  (25) 

X z 11  a 


h/2 


^y'^z  " ^12  a 


(26) 


where  the  effective  diffusivities,  d^^  and  d^^  ^^re  defined  by  the  relationship 

h/2 


= i / 


(d,I,dl2)d. 


(27) 


-h/2 


A similar  problem  with  the  finite  dimension  in  the  y direction  yields 


22 


i / 

-h/2 


h/2 


d_-d 
22  z 


(28) 
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SECTION  V 

SOLUTION  FOR  A LAMINATED  PLATE 


i 

I 

i 


In  most  practical  applications,  fiber  reinforced  composite  laminates 
are  constructed  of  equal  numbers  of  plies  at  ±6.  A cursory  examination  of 
equation  (4)  shows  that 

dj2(-W)  = - (30) 

Equation  (30)  in  conjunction  with  equation  (29)  for  an  equal  number  of  ±0  plies 
yields 

dj2  = 0 (31) 

Thus  for  most  practical  laminates,  the  governing  equation  for  diffusion  is  of 
the  form 

m . = d,,m  + d__m  +d,,m  (32) 

,t  11  ,xx  22  ,yy  33  , zz  ' ’ 

Consider  a rectangular  laminated  plate  containing  an  initial  uniform 
moisture  content,  m^,  with  the  origin  of  the  coordinate  system  at  the  center 
of  the  plate,  and  subjected  to  a constant  temperature  and  humidity  condition, 
as  described  by  the  following  initial  and  boundary  conditions 

m = m^  , t<0  (33) 

m(d:a/2,y,  z)  = m(x,:feb/2,  z)  (34) 

= m(x,  y,  dh/2)  = m 

A solution  to  equation  (32)  which  satisfies  the  conditions  of  equations  (33)  and 
(34),  is  of  the  form  (Reference  4) 

m = m +(m  - m )g  (35) 

o e o 

where 


g = 1 - mj(Xj,t)m2(x2,t)m2(x2,t) 


I.  ,,.^11  


and  m.(x.,t)  are  solutions  to  the  one  dimensional  equation 


m = d..m  , i = 1,2,  3 

i,t  ii  »x.3c. 


(37) 


subjected  to  the  conditions 


m.(x,  , t)  = 1 t < 1 


m.(±a^/2,t)  = 0 , t >0 


(38) 

(39) 


where  x^,  x^,  x^,  and  a^,  a^,  a^  correspond  to  x,  y,  z,  and  a,  b,  h, 
respectively. 


Integration  of  equation  (34)  over  the  plate  volume  yields 


M = M + (M  - M )G 
o e o 


(40) 


where 


G = 1 - Mj(t)M2(t)M3(t\ 


(41) 


and 


a,/2 


M.(t)  = 


J m.(x.,t)dx. 


(42) 


-aj/2 


1 . TRIGONOMETRICAL- SERIES  SOLUTION 


A solution  to  equation  (37)  which  satisfies  the  conditions  of  equations 
(38)  and  (39)  can  be  obtained  by  classical  separation  of  variables  with  the 
result  (Reference  1) 


m. 

1 


4 ^ 


= ? Z 


rrx. 


k=0 


(2k+l) 


cos(2k+l) 


(43) 


x exp  [-(2k+l)^  if^t*  ] 


J 


where 


* 

t. 

1 


d..t 

n 

2 

a. 

1 


Integration  of  equation  (43)  over  the  plate  volume  yields 


M = ^ — z [-(2k+l)^fr^t*] 

” k=0  (2k+l) 


2.  LAPLACE  TRANSFORM  SOLUTION 


(44) 


(45) 


Short  time  solutions  for  either  total  weight  gain  or  moisture  distribu- 
tion are  often  of  practical  interest.  Because  of  the  slow  convergence  of  equa- 
tions (43)  and  (45),  an  alternate  solution  to  equations  (37),  (38),  and  (39)  which 
provides  rapid  convergence  for  short  times  is  desirable.  Such  a solution  can 
be  obtained  from  the  Laplace  transform  method  with  the  result  (See  Appendix) 


Integration  of  equation  (46)  yields 


M. 

1 


1 +2Vr  £ (-1)*" 

k=l 


X ierfc 


k 


(47) 


3.  APPROXIMATE  SOLUTION 

It  has  been  noted  by  Weitsman  (Reference  3)  that  for  t.*<  0.01,  equation 
(46)  can  be  accurately  represented  by  retaining  only  the  first  term.  For 


10 


t.*  < 1,  the  first  three  terms  are  sufficient  for  an  accurate  answer.  Thus, 
the  three-dimensional  solution  can  be  accurately  approximated  by  the  equations 


g = 1 - erf 


i-2x/a 


* (izl^ 


t*<  0.01 


g = 1-erf 


l-2x/a 


l-2y/b 


1 - z (-»' 

k=l 


+ erfc 


2k- 1 +2z/h 


2k-l-2z/h 


0.01  < t*<  1 


where 


Equations  (48)  and  (49)  utilize  the  identity 


1 - erfc  u 


erf  u 


Equations  (49),  (50),  and  (51)  are  valid  for  R , R < 0.02. 

1 M 

For  t.  < 0.07,  equation  (47)  can  be  adequately  represented  by  the  first 

term  of  the  series  (Reference  1).  Similarly,  a cursory  examination  of 

% 

equation  (45)  reveals  that  for  t^  >0.07,  the  first  term  is  adequate.  Thus, 
for  Rj,  R^  < 0. 1, 


t^<  0.07 


o = .-i 


1 - 4 


X exp  ( - »r  t*)  , 0 . 07  < t*  < 1 

G = 1 , t*  > 1 


4.  SHORT  TIME  SOLUTION 


Equation  (53)  in  expanded  form  yields  a third  order  polynomial  inv/tT” 
For  R , R_  <0.1  and  t*<0.01,  only  the  first  order  term  need  be  retained 
with  the  result 


G = 


where 


and  d is  an  effective  diffusivity  defined  by  the  relationship 

d = dj^d  +\/r^ 

Thus  a plot  of  G versus  VT~  yields  a straight  line  over  the  initial  portion  of 
the  curve.  The  slope  of  the  linear  region  is  directly  related  to  d,  i.e. 


d = 


Equation  (59)  provides  a relationship  for  experimentally  determining  d. 


5. 


ONE  DIMENSIONAL  APPROXIMATION 


For  the  case  where  R^,  R^  < 0.002,  equation  (56)  can  be  used  in  con- 
junction with  a one  dimensional  solution  for  total  weight  gain  to  predict  an 
approximate  three  dimensional  total  weight  with  the  result  that  equations 
(52),  (53),  and  (54)  become 


fr* 

.It  —* 

= 4^—  . t <0. 


07 


8 2^^ 

G = 1 - —T  exp  (-  IT  t ), 


0.07<  t*<  1 


— * 

t > 1 


(60) 

(61) 

(62) 


This  approximation  was  first  suggested  by  Shen  and  Springer  (Reference  6). 

Note  that  for  t = 0.07,  equations  (60)  and  (61)  yield  a value  of 

/— ♦ 

G = 0.6.  Thus,  the  total  weight  gain  will  be  a linear  function  o£\  t up  to 
approximately  60%  of  equilibrium. 


SECTION  VI 
NUMERICAL  RESULTS 


Consider  a laminate  constructed  from  plies  of  the  same  unidirectional 
material  having  symmetric  fiber  packing.  As  a result  d^^  = d^.  For  sim- 
plicity, consider  laminate  stacking  geometries  which  yield  = ^2Z*  ^ 
cursory  examination  of  the  transformations  given  in  equation  (4)  shows  the 
inplane  diffusion  coefficients  to  be  independent  of  0 for  this  case.  Thus,  the 
laminate  is  qua  si- isotropic  with  respect  to  the  inplane  diffusivities . Integra- 
tion of  equation  (4)  over  all  possible  values  of  0 yields  the  following  results 
for  quasi- isotropic  diffusion 

i i L T , / 


Using  equation  (7)  in  conjunction  with  equation  (63)  yields 


1-  - Vj/2 


The  volume  fraction  is  chosen  as  0.6  for  the  present  report.  Actual  volume 
fractions  for  graphite /epoxy  laminates  vary  from  about  0.6  to  0.65.  Equation 
(64)  now  yields  the  result 


= 2.09 


In  Figure  2 the  moisture  profile  for  two  dimensional  diffusion  is  shown 

* 

across  the  centerline  of  a quasi-isotropic  laminate  for  two  values  of  t . 
Numerical  results  for  total  weight  gain,  along  with  experimental  data,  is 
shown  in  Figure  3 for  a unidirectional  graphite /epoxy  composite  and  in 
Figures  4 and  5 for  graphite /epoxy  laminates  having  different  thicknesses. 
Equations  (7),  (51),  (54)  and  (65)  were  used  in  conjunction  with  specimen 
dimensions  to  generate  the  theoretical  curves. 


SECTION  vn 
EXPERIMENTAL  DATA 


In  order  to  assess  the  accuracy  of  Fick's  law,  composite  plates  were 
fabricated  from  Hercule's  AS/ 3501 -5  graphite /epoxy  pre-preg  system,  cut 
into  small  specimens,  and  exposed  to  equilibrium  under  various  temperature 
and  humidity  conditions.  One  panel  was  a four  ply  unidirectional  from  which 
2 in.  (5cm)  x 2 in.  (5cm)  specimens  were  cut.  Two  laminate  panels  were  also 
fabricated.  One  panel  was  a four  ply  (0,  90)  construction  from  which  2 in. 
(5cm)  X 2 in.  (5cm)  specimens  were  cut.  The  second  panel  was  constructed 
of  twenty  plies  with  the  stacking  geometry  (0,  90)^^  from  which  1 in.  (2.54cm) 

X 1 in.  (2.54cm)  specimens  were  cut.  Humidity  exposures  included  underwater, 
UW,  95%,  and  75%.  Temperatures  included  100°F  (38°C),  120°F  (49°C), 

140°F  (60°C),  160°F(71°C).  Exposure  procedures  were  in  accordance  with 
those  described  in  Reference  5. 

For  the  absorption  data,  percent  weight  gain  was  recorded  for  various 
exposure  times.  A value  of  M^  was  determined  by  denoting  an  apparent 
equilibrium  moisture  content  from  the  data.  The  weight  gain  values  were  an 
average  of  at  least  three  specimens. 

Values  of  M/M^  for  the  unidirectional  c mposite  specimens  were  tabu- 
lated as  a function  ofVt” for  each  exposure  temperature.  A linear  regression 
analysis  was  used  to  fit  the  data  over  the  range  0^  M/M^<.  0.6,  since  R^, 

R_  < 0.002  for  the  dimensions  of  the  unidirectional  specimens.  The  slope  of  this 

to 

line  was  then  used  in  conjunction  with  equation  (59)  to  determine  d.  For 

absorption  from  an  initially  dry  condition,  = 0,  equation  (40)  reduces  to 

G = M/M  . 

e 

Equation  (58)  in  conjunction  with  equations  (7)  and  (51)  yields 
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SECTION  VIII 
DISCUSSION 


The  importance  of  accelerated  environmental  conditioning  techniques  is 
illustrated  in  Figure  2,  In  particular,  the  severe  changes  in  moisture  gradient 
near  the  edge  with  respect  to  t*  are  independent  of  material  diffusion  properties, 
but  the  rate  of  change  of  these  gradients  depends  on  the  diffusion  rates.  Thus, 
high  temperature  moisture  conditioning  which  produces  large  diffusion  coeffi- 
cients may  induce  edge  cracks  because  of  the  rapid  change  in  the  moisture 
gradient  at  the  edge. 

A comparison  of  the  data  with  Pick's  law  is  illustrated  in  Figures  3-5. 

A cursory  examination  of  equations  (53)  and  (54)  reveals  that  G is  determined 
directly  by  the  value  of  t*  if  and  remain  constant.  Thus,  a plot  of 
M/M^  versus  \j t*  allows  data  from  different  relative  humidities  and  tempera- 
tures to  be  plotted  on  one  master  plot. 

Good  correlation  between  data  and  theory  is  obtained  for  the  unidirectional 
composite  and  thick  bidirectional  laminate  while  the  correlation  for  the  thin 

bidirectional  laminate  is  only  fair.  A linear  regression  analysis  of  the  data  in 

• • ^ 

Figure  4 yields  a good  fit  to  the  theory  for  small  values  of  t with  the  resulting 

values  of  d,^,  being  considerably  larger  than  those  listed  in  Table  1 for  the  uni- 
directional composite.  This  procedure,  however,  leads  to  a poor  fit  for  large 
* 

values  of  t . This  anomaly  has  been  discussed  in  Reference  5 with  respect  to 
a possible  cracking  phenomenon  or  stress  dependent  diffusion  process  associated 
with  large  residual  stresses  which  are  present  in  bidirectional  laminates. 

The  fact  that  the  same  anomaly  does  not  appear  to  occur  in  the  thick  laminate 
makes  the  precise  mechanism  more  difficult  to  determine. 
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Experimentally  determined  values  of  d^  are  shown  in  Table  1 as  a function  of 
temperature. 


Table  1 

EXPERIMENTALLY  MEASURED  DIFFUSIVITIES 


T 

d^(IN^/HR) 

d^(MM^/SEC) 

100°F(  38°C) 

5.65  X lO'^ 

1.01  X 10"^ 

120°F(49°C) 

8.21  X lo'^ 

1.47  X lO"^ 

160°F(71°C) 

2.02  X 10‘^ 

3.62  X lO'^ 

180°F(82°C) 

3.91  X lO"^ 

7.00  X 10“^ 

These  values  of  d were  used  in  conjunction  with  experimentally 

* 

measured  exposure  times  to  establish  values  of  t corresponding  to  the 
measured  weight  gains  for  the  bidirectional  laminates. 


SECTION  DC 

CONCLUDING  RElvLARKS 


The  mathematics  of  diffusion  in  a laminated  composite  can  be  described  I ; 

J I ' 

by  simplified  closed  form  expressions  for  the  case  of  constant  temperature  **  i 

and  constant  humidity  exposure.  These  models  assume  the  diffusion  process 

to  be  independent  of  concentration.  Comparison  to  data  gives  some  confidence 

in  the  mathematical  models.  Anomalies  associated  with  thin  bidirectional  : 

laminates  need  to  be  resolved.  Stress/strain  effects  on  the  diffusion  process, 

time  dependent  material  property  changes,  and  cracks  due  to  swelling  are  items  ' 

which  need  further  investigation.  The  micromechanics  of  diffusion,  equations 

(5)-(7),  do  not  recognize  the  interface  as  a possible  source  of  diffusion.  A i 

more  precise  definition  of  the  mechanisms  associated  with  diffusion  behavior  | 

are  necessary  before  modified  or  advanced  mathematical  models  of  the  | 

diffusion  process  can  be  developed  at  either  the  micromechanical  or  macro-  ] 

mechanical  level.  j 
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